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UNSTABLE  CAPILLARY  WAVES  ON  SUEFACE  OF  SEPAEATION 

OF  WO  VISCOUS  FLUIDS* 

By  V.  A.  Borodin  and  Y.  F.  Dityakin 


The  study  of  the  breakup  of  a  liquid  jet  moving  in  another  medium, 
for  example,  a  jet  of  fuel  from  a  nozzle,  shows  that  for  sufficiently 
large  outflow  velocities  the  jet  breaks  up  into  a  certain  number  of 
drops  of  different  diameters.  At  still  larger  outflow  velocities,  the 
continuous  part  of  the  jet  practically  vanishes  and  the  jet  immediately 
breaks  up  at  the  nozzle  into  a  large  number  of  droplets  of  varying 
diameters  (the  case  of  "atomization") .  The  breakup  mechanism  in  this 
case  has  a  very  complicated  character  and  is  quite  irregular,  with  the 
droplets  near  the  nozzle  forming  a  divergent  cone. 

Eayleigh  (reference  l)  was  the  first  to  make  a  theoretical  study 
of  the  jet  and  to  establish  the  possibility  of  droplet  formation.  The  ■ 
disturbance  of  a  jet  of  an  ideal  fluid  flowing  into  a  vacuum  and  having 
a  wave  length  4.4  times  as  large  as  the  diameter  of  the  jet  is  shown 
to  grow  more  rapidly  than  other  disturbances;  eventually,  the  jet 
breaks  up  into  droplets  of  the  same  diameter.  Eayleigh  succeeded  in 
determining  theoretically  the  drop  diameter,  the  value  of  which  agrees 
well  with  teats  on  jets  issuing  with  very  small  velocities.  Later, 
the  viscosity  of  the  jet  was  also  taken  into  consideration.  The 
viscosity  is  found  to  decrease  the  rate  of  amplitude  increase  of  the 
disturbances  but  the  ratio  of  the  optimal  length  of  the  wave  to  the 
diameter  of  the  jet  remains  unchanged. 

Other  authors  that  studied  the  conditions  of  the  axial-symmetrical 
breakup  of  a  jet  of  a  viscous  liquid  found  that  the  ratio  of  the  optimal 
wave  length  to  the  jet  diameter  was  somewhat  greater  than  that  computed 
by  Eayleigh. 

In  addition  to  the  viscosity,  Tomotika  (reference  2)  took  into 
account  the  density  and  viscosity  of  the  medium  surrounding  the  jet 
and  obtained  good  agreement  with  tests  on  jets  issuing  with  very  small 
velocities  for  which  droplets  of  the  same  diameter  are  formed. 


*"Neustoichivye  Kapilliarnye  Volny  na  Poverkhnosti  Eazdela  Dvukh 
Vyazkikh  Zhidkostei."  Prikladnaya  Matematika  i  Mekhanika.   Vol.  XIII, 
no.  3,  1949,  pp.  267-276. 
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Neither  of  the  aforementioned  theories  of  the  breakup  of  a  liquid 
jet  provided  a  basis  for  the  phenomenon  for  the  case  of  breakup  into 
droplets  of  different  diameter,  a  fact  that  is  explained  by  the 
idealized  conditions  of  the  problem.  This  idealization  consisted  either 
in  neglecting  the  viscosity  of  the  jet,  the  density,  and  viscosity  of 
the  surrounding  medium,  or  the  inertial  forces.  Such  simplifications 
were  assumed  in  view  of  the  complicated  mathematical  equation  (generally 
transcendental)  that  determines  the  relation  between  the  wavelength 
and  the  increment  of  the  vibration  amplitude. 

In  the  present  paper,  an  attempt  is  made  to  provide  a  mathematical 
basis  for  the  possibility  of  the  appearance  of  droplets  of  different 
diameters  as  a  result  of  the  jet  breakup  on  the  basis  of  the  considera- 
tion of  unstable  capillary  waves  on  the  surface  of  separation  of  two 
viscous  liquids. 

For  simplification  of  the  solution  of  the  problem,  particularly 
for  obtaining  the  algebraic  characteristic  of  the  equation,  the 
lengths  of  the  capillary  waves  on  the  surface  of  the  liquid  jet  are 
assumed  to  be  so  small  in  comparison  with  the  jet  radius  that  the  jet 
may  be  considered  infinitely  large;  study  of  the  stability  of  the 
plane  surface  of  separation  of  two  infinitely  extending  viscous  fluids 
can  thus  be  made.  This  assumption  represents  a  considerable  degree  of 
idealization  but  nevertheless  permits  a  qualitative  explanation  of 
not  one  but  several  unstable  capillary  waves  that,  in  passing  through 
the  jet,  lead  to  the  formation  of  droplets  of  differing  diameters. 

The  existence  of  several  unstable  capillary  waves  is  demonstrated 
that  can  lead  to  the  breakaway  of  several  infinitely  long  strings  of 
different  dimensions  from  the  partition  surface.  The  problem  investi- 
gated gives  a  rough  approximation  of  the  disintegration  pattern  of  a 
liquid  jet  in  another  medium  and  does  not  pretend  to  explain  the  com- 
plicated mechanism  of  the  limiting  form  of  the  disintegration  of  a 
jet,  namely,  atomization.  Nevertheless,  one  of  the  peculiarities  of 
atomization,  the  appearance  of  a  dimension  spectrum  of  the  droplets, 
begins  to  appear  even  for  the  given  idealized  consideration  of  the 
stability  of  the  partition  surface. 

1.  Equations  of  small  waves  and  their  solution.  -  A  plane  surface 
of  separation  of  two  infinitely  extending  viscous  fluids  (fig.  l)  is 
considered.  The  viscosity  and  density  of  the  lower  fluid  are  denoted 
by  i-i-L  and   p1;  respectively,  and  of  the  upper  fluid  by  u2  and  p2* 
The  lower  fluid  is  a.ssumed  to  move  with  the  velocity  V-j.  an&   the  upper 
fluid  with  the  velocity  V2,  the  direction  of  motion  being  the  same  and 
the  velocities  independent  of  y. 
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A  study  of  the  character  of  the  equilibrium  of  the  surface  of  sep- 
aration under  the  action  of  the  viscous  forces  and  the  forces  of  surface 
tension  that  impart  to  both  liquids  small  disturbances  parallel  to  the 
x-axis  is  presented.  The  fluids  shall  be  considered  incompressible  and 
weightless  and  shall  cause  certain  disturbances  to  the  components  of  the 
motion. 


vx  =  V  +  ^x 


V  =  v 

y   y 


P  =  p  +  P^ 


It  is  further  assumed  that  the  velocities  of  the  imposed  disturb- 
ances and  their  derivatives  up  to  the  third  inclusive  are  small  and  that 
the  magnitudes  of  the  second- order  smallness  may  be  neglected. 


From  the  Navier-Stckes  equations,  the  following  equations  of  the 
imposed  disturbances  are  obtained: 


dvx 
~5T 


+  V 


avx 


1  dpj*; 
P  ox 


+  uAv^ 


oV      ovy 
ot      ox 


P  oy      y 


(1.1) 


where   u  =  |j/p  is  the  kinematic  viscosity. 
The  equation  of  continuity  is 


Sv-^ 


0 


oX    By 
By  introducing  the  stream  function  of  the  disturbance 


(1.2) 


By 


ox 


(1.3) 


and  by  eliminating  the  pressure  p*  from  equations  (l.l),  the  idealized 
equation  is  thus  obtained  in  the  Helmholtz  form 


oX 


at 


(1.4) 
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Let  the  stream  function  of  the  imposed  disturbance  be  a  periodic 
function  of  x  and  of  the  time  t: 


„  =  f(y)ei(-et)        £..*)  (1.5) 


where  a  is  the  propagated  circular  frequency  of  the  vibrations  (the 
wave  number),  A.   is  the  wavelength  of  the  imposed  disturbance, 
p  =  Pr  +  ip^  is  the  complex  frequency  of  vibrations  in  time,   P 

is  the  real  frequencj^  of  vibration  in  time,  and  p^  is  the  increment 

of  the  growth  of  vibration  or  the  decrement  of  damping. 

The  character  of  the  wave  motion  on  the  surface  of  separation 
after  the  imparting  of  disturbances  to  both  surfaces  will  thus  depend 
on  the  sign  of  the  imaginary  part  of  the  frequency  0j..  If  Pi   is 

positive,  there  will  be  an  increase  in  the  wave  amplitude  with  time; 
if  Pi  is  negative,  there  will  be  a  damping  of  the  wave  amplitude; 
finally,  if  Pr  =  0,  there  will  be  an  aperiodic  increase  (0^  >  0) 

or  a  damping  (Pj  <  0)  of  the  wave  amplitude.  By  substituting  expres- 
sion (1.5)  in  equation  (1.4),  the  following  equation  is  obtained: 

UfIV  -  (2co2U  -  10)  f "  -  (ipa2  -  pa4)  f  -  iVa  (f "  -  a2f)  =  0     (1.6) 

The  problem  of  the  characteristic  values  of  a  homogeneous  system 
of  equations  of  the  fourth  order  will  be  considered. 

By  setting  f "  -  ccrf  =  tp,  a  system  of  equations  of  the  second 
order  is  obtained. 

<p«+  (l£~S.  a2V  =  0     f "  -  a2f  =  q>         (1.7) 

Hereinafter,  the  following  notations  are  introduced: 


J 


P  -  vVl  /  0  -  V2a 


a2  =  m     A/i  — a^  =  m9  (1.8) 

1  v     v2 


The  solution  of  the  first  of  equations  (1.7)  has  the  form 

qp  =  (^e1™!^  +  Cge-^iy  (1.9) 

By  substituting  expression  (1.9)  in  the  second  of  equations  (1.7), 
a  non-homogeneous  equation  is  obtained  for  which  the  solution  is 
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-F  = 


.imiy 


w-,2  +  a2 


C'l 


■lnny 


2  2      "2 

m-,      +  or 


Co   +  eay  C,   +  e"ay  C, 


(1.10) 


The   stream  function  for  the   lower  and  upper  liquids  according  to 
equation    (1.5)  will  be 


i(ox-pt) 


imiy 


■imiy 


m-,      +  cc' 


2      1 


m-,      +  a' 


2   C2  +  e"*C3+e-«*C 


(1.11) 


V  2=  e 


i(cac-Pt) 


,in»2y 


-im2y 


mo     +  a 


Bio     +  ar 


Cc  +  eay  C-7   +  e 


-ay 


(1.12) 


The  arbitrary  constants  C^  must  be  determined  from  the  conditions 
on  the  surface  of  separation  and  at  infinity. 

2.  Boundary  conditions.  -  The  boundary  conditions  of  the  problem 
will  be  as  follows: 

1.  At  infinity  (y  =  +<*>),   finite  solutions  must  be  maintained  for 
^-j_  and  '4<  2'  Hence,  the  arbitrary  constants  of  the  terms  with  positive 

exponents  for  ty  n   and  with  negative  exponents  for  \|/  o.     must  be  equated 

to  zero:  C-,  =  C3  =  Cg  =  Cg  =  0.  Thus,  equations  (1.10)  and  (l.ll)  will 

have  the  form. 


¥]_  =  e 


i  (ax-(3t] 


^2  =  e 


=   i(ax-pt)  f _ 


—2 2  C2  +  e  ^  Ca 

m-^  +  ar 


(2.1) 


m2  +  ar 
On  the  surface  of  separation,  there  must  be  no  slip,  that  is, 


(V.Tl)y=0"  My=0 


or 
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o^ 


St; 


ey/y=0       \oy/y=0 


2 

ox; 


y=o 


(2.2) 


3.  The  tangential  stresses  on  the  surface  of  separation  are 
continuous 


^AViV=0  =  ^AV2)y=0 


(2.3) 


4.  The  difference  between  the  normal  stresses  p  -j_  and  p  g 

on  the  surface  of  separation  is  equal  to  the  pressure  brought  about  by 
the  surface  tension;  that  is, 


(- 


dv. 


b2h 


"?yl  "  %2  =  [~   ?!  +  2^l  -£r  I    -    I  "  p2  +  2^2  ^&  J  =  -  T  -__2     (2.4) 


ay 


oy 


ox 


where  T  is  the  capillary  constant  of  one  liquid  relative  to  the  other 
and  h  is  the  rise  in  the  surface  of  separation  at  the  point  x. 

By  using  equations  (2.1),  the  boundary  conditions  (2.2)  are  obtained 
in  the  form 


im]_ 


imc 


2  ,  „2  C2  ~  ^4  +  „  2  .  „2  u5  "  ^7 


aC7  =  0 


ml 


m-p  +  a 


^2  ^5 

2 2  +  C4  +  — 5    2  _  C7  =  ° 

m-^  +  a.^        nu^  +  ar 

Similarly,  the  boundary  condition  (2.3)  is  obtained  in  the  form 


(2.5) 


^2 


C2 


m-i   +  a 


+  CA\a?   +  — ~ -7T   C2  +  a2C4 

m-|  6   +   ar 


mo  +  a 


+  C7  \a?   + 


m9  9 

— S Cq  +  a^C7 

2  ,  „2   5      7 


(2.6) 


mg   +  cc 


NACA  TM  1281 


The  pressures  p^  and  P2  are  computed  from,  equations  (l.l) 
and  (2.1).   Thus 


Pi  =  Pie 


i(ox-Bt) 


ae 

mn 


I f  _± —  _  — i_  +  iR  +  m-,  V-,  \  C2  + 

^  +  a2  V  ml      ml    ml       J 

(p  +  iU  a2  -  V  a  -  il)a2e-ay  C^ 


(2.7) 


2  =  p2e 


Pq  = 


(CQC-Pt) 


^im2J  /  uo<x    iVoa,   ,•  r 


m22  +  a2  V  m2 


nu    mo 


The  rise  of  a  point  on  the  surface  is  a  periodic  function  of  x 
and  t. 

h=Hei(ax-Pt)  (2>Q) 

where  H  is  the  maximal  rise  of  a  point  on  the  surface  of  separation. 
The  velocity  of  the  raised  point  on  the  surface  of  separation  is 


(vyiL=o  = 


^i 


m 


+  V 


y=u    Vdx/y=o     ot 


1 


Sh 
ox 


(2.9) 


After  differentiating  expressions  (2.1)  and  (2.8)  and  by  substi- 
tuting in  expression  (2.9),  the  following  equation  is  obtained: 


H  = 


a 


aV,  -  3 


C„  - 


(2.10) 


m-i   +  a 


By  substituting  equation  (2.10)  in  (2.9)  and  by  differentiating 
equation  (2.9), 


cr- 


o2h 

ix2  "   "i  -   U2   , 


2-^  -  C/U1^-^) 


(2.11) 
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By  computing  the  derivatives  dv  ,/dy  and  6v  2/dy  and  substi- 
tuting them  simultaneously  with  expressions  (2.7)  and  (2.11)  in  (2.4), 
the  following  boundary  condition  is  obtained: 


a 


m-i      +  of1 


'^l0-    +  iplVla  "    i'3pl  Ta' 
+  m1|i1  + 


m 


1 


aV1  -   B 


Co   - 


p-j_P   -   P^a  -   i2u1a2   + 


aYn    -   B, 


C„    + 


a 


2  2 

m.0     +  ar 


\±2a    ~  ip2^2a  +   ^p2 


m- 


3m2u2 


C5   +    ^p2|3   "    p2V2a  +   i2M2a  )    c7   =  ° 


The  following  nondimensional  parameters  are  then  introduced 


(2.12) 


Z  = 


U1a 

Rl 

Vl 
V 

K2 

V2 

A 

Ul 
~U2 

N  = 

Tpx 

2 

\i1  a 

K 

.  iii 

(2.13) 
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where  c  =  p/cc.  is  the  complex  wave  velocity. 

Equations  (1.8)  can  then  he  represented  in  the  forms 


mi  =  claJ±(Z  -   E1)  -  1 


=  a/^i(ZA  -  E2)  -  1 


Equations  (2.5),  (2.6),  and  (2.12)  are  represented  in  nondimen- 
sional  parameters.  The  following  notations  are  first  introduced: 


&1*  =  | 


1  Z-E-, 


N 


A/i  (Z-K-,  )    -  1       Z~E1 


=  ha 


l"! 


■:     =  a\  (  Z-B1-21   +  -2-\  =  cc2^ 


'1. 


R1"Z> 


2(i   -   2E2   +  2AZ)  „ 

C-i*    =    Uo    ; ■  =    UoCrCn 


a3*=--^      j 


a2  Z_E1  a2 


dx*  =  u2a2    (ZA-E2   +  2i)   =  u2a2d1 


1        i 


a 


2   ZA-E 


2        a? 


,*  =  i  a/i(Z-R,  )    -   1  =  — 
2  a    V     v        1'  a 


c?*  =  - 


i  aJUz^Bz)  -  i      c2 


ZA-S. 


a 


>  (2.14) 


where     a-^,   a2,   aj,  b-j_,   c-j_,    c2,    c-j,    and  d^     are  likewise  nondimen- 
sional  magnitudes. 
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The  following  system  of  equations  is  then  obtained  for  the  con- 
stants C2,  C4,  C3,  and  C^: 


a1*C2   -  b1*C4  +  c-j*C5  +  d1*C7  =  0 


-i 


a2*C2  " 

aC4  +  Cg*C5  -  a 

c7  =  0 

ag*C2  + 

C4  +  c3*C5  - 

c7  =  0 

K  C2 


=  0 


(2.15) 


This  system  of  homogeneous  equations  has  solutions  different 
from  zero  if  its  determinant  is  equal  to  zero.  By  setting  up  the 
determinant  and  expanding 

2K(a1  +  c1)  +  (dx  -  Kb1)(a2  +  Kc2)  +  (Kb]_  +  dx)  (Kc3  -  a3)  =  0 

By  solving  this  equation  for  Z,  the  following  wave  equation  of 
the  18th  degree  with  complex  coefficients  is  obtained: 


r18Z18  +  (rin   +  is,,)  Z1' 


,17 


17 


17' 


..*  +  {r1   +  is1)  Z  +  (rQ  +  isQ)  =  0 

(2.16) 


The  real  and  imaginary  parts  of  the  coefficients  depend  on  the 
five  nondimensional  parameters:  R-j_,  R2,  A,  N,  and  K. 


5.  Investigation  of  roots  of  characteristic  equation.  -  The 
increase  in  oscillation,  that  is,  the  loss  of  stability  of  the  sur- 
face of  separation,  arises  from  those  waves  for  which  the  imaginary 
part  of  the  frequency  is  positive  (,Bj  >  0).  Hence,  the  investiga- 
tion of  the  roots  of  equation  (2.16)  should  determine  those  ranges 
of  the  parameter  N  or  the  wave  number  a  in  which  the  complex 
roots  of  the  equation  lie  in  the  upper  half -plane. 

By  the  Rayleigh  hypothesis,  the  further  development  of  an 
unstable  deformation,  that  is,  the  form  and  dimensions  of  the  parts 
breaking  away,  is  determined  by  the  critical  (or  optimal)  disturb- 
ances.  The  critical  disturbances  may  be  defined  as  those  that 
develop  more  rapidly  than  the  others  or  that  correspond  to  the 
maximum  increment  of  the  growth  $±.     This  principle  of  deter- 
mining the  character  of  the  unstable  deformations  by  the  character 
of  the  maximum  unstable  disturbance  has  been  experimentally 
confirmed  by  a  number  of  investigators  (reference  3). 
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In  the  case  considered,  the  growth  in  the  amplitudes  of  the 
oscillations  will  lead  to  breakaway  of  infinitely  long  strings 
from  the  surface  of  separation,  similar  to  the  formation  and 
breakaway  of  wave  crests.   The  separation  will  take  place  for  such 
values  of  a  or  wavelengths  k     for  which  p^  has  the  maximal 
value . 

If  a  spectrum  of  small-period  disturbances  that  can  be  developed 
into  a  Fourier  series  can  be  assumed  to  be  imposed  on  both  liquids, 
the  harmonics  with  the  wavelengths  equal  to  the  wavelengths  of 
the  maximal  unstable  disturbances  bring  about  a  separation  of 
infinitely  long  strings  from,  the  partition  surface.   Because  the 
characteristic  dimension  (for  example,  the  diameter  of  the  transverse 
string)  is  connected  with  the  length  of  maximal  unstable  disturbance, 
strings  of  different  dimensions  will  break  away  from  the  surface  of 
separation.   In  figure  ?, ,   the  scheme  of  formation  of  such  strings 
for  three  successive  instants  of  time  is  shown. 

Investigation  of  the  roots  of  the  simplest  particular  case  of 
equation  (2.16)  is  presented. 

Let  both  fluids  be  stationary  and  their  kinetic  viscosities 
the  same.   In  this  case,  V-j_  =  Vg  =  0,  u^  =  v^}   m.-.  =  nig,  A  =  1, 

R-j_  =  Kg  =  0,  and  equation  (2.16)  goes  over  into  an  equation  of  the 

8th  degree  whose  coefficients  depend  only  on  the  two  parameters 
K  and  N: 

AQZ8  +    (A]    +   iB1)Z7   +    (A2   +  iB2)Z6  +   (A     +  iB3)Z5  +    (A     +  iB4)Z4  + 

(A5  +  iB5)Z3  +    (A6  +  iB6)Z2   +   (A7   +  iB?)Z     +  A8   =  0  (3.1) 
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where 


AQ  =  -  (1  -  K)2 

A2  =  2E  (K  -  1)  N  -  E4  +  2E3  -  4K2  +  6K  +  13 

A-,  =  -  2K  (K  -  I)2 

A3  =  4E2  (K  -  1)  W  -  2K  (3K2  +  13) 

A7  =  2E3K2 

A4  =  -  E2E2  +2  (E4  -  E3  +  5E2  +  5K)  W  -  12K3  +  26K2  -  10K  -  9 
A5  =  -  2E3H2  +  12E3N  -  8K3  -  8K2 

A8  =  -  K2  (1  +  2K)  N2 

=  (1  -  K2)  E2N2  +  (4E4  -  10E3  -  4K2  -  6K)  N  -  SE3  +  12E2  -  8E  +  4 

S-l  =  2  (E2  +  2E  -  3) 


B 


o 


3E4  -  14E3  +  13E2  +  18E  +  13  -  SEN 

B2  =  4E  (E  -  I)2 

B4  =  3K3  -  4E2  -  20E  -  2E3K 

B6  =  4E2  (1  -  E)  IT 

B5    =   -   2E%2   +   [2E    (1   +  K)    (1  +  E  -  E2)    +  SE3   +  4E2   +   4E]n  + 

4(E   -    1   -  E2)    (l   +  E    -  E2)    +  4E4    -   20E3    -    4(K  -    1    -  E2)2 

B7   =   [E2    (1   +  E)2   -   2E4]n2   +    [4E4  +  4E    (l  +  E)    (E  -   1   -  E2)] 


(3.2) 
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The  charact eristic  equation  (3.1)  is  a  polynomial  vhose 
coefficients  depend  nonlinearly  on  the  two  parameters  K   and  N. 
Each  pair  of  values  of  the  parameters  K  and  K  or  each  point 
of  the  plane  EN  correspond  to  the  completely  defined  polynomial 
(3.1),  that  is,  completely  determined  values  of  the  eight  roots 
of  the  polynomial.  In  the  plane  Etv,  it  is  evidently  possible  to 
find  a  curve,  each  point  of  which  corresponds  to  the  polynomial 
(3.1),  that  has  at  least  one  root  located  on  the  real  axis  so  that 
only  in  crossing  this  curve  is  a  crossing  of  the  roots  through  the 
real  a,xis  possible.   This  curve  "breaks  up  the  plane  EN  into 
regions,  the  points  of  which  each  correspond  to  polynomials  (3.1), 
that  have  the  same  number  of  roots  with  positive  imaginary  part. 

These  curves  are  constructed  by  making  use  of  the  method  of 
Y.  I.  Neimark  (reference  4)  that  permits  a  breakup  of  the  plane 
of  the  parameters  for  the  roots  of  the  polynomial  lying  in  the  left 
or  right  half -plane. 

The  substitution  Z  =  -it,     is  made.  The  upper  half -plane  of 
the  roots  of  equation  (o.l)  is  transformed  into  the  left  half -plane 
of  the  roots  of  the  equation 

-  AQg   +   i(Ax   +  iB^C7   +    (A2   +  iB2)t6   -   i(A3   +  iB3)$5  -    (A4   +   iB4)£4   + 

i(A5   +  iB5)£3     +    (A6   +   iB6)C2      -    (A7   +   iB7)C     +  A8   =   0 

(3.3) 

By  substituting   £  =  i£Al   in  the  preceding  equation  and  multi- 
plying the  result  by  tjS,  equation  (3.3)  is  reduced  to  the  form 

F(£,ti)  +  iG(£,n)  =  0  .         (3.4) 

where 

F(t,n)    =  A0e8   +  Axe7Tl   +  A2|6n2   +  A3£V  +  A4|4n4  +  A5|3n5   + 

A6e2T!6   +  A-^7   +  A6ti8  (3.5) 

G(U)    -  ^fn    +  B2e6ri2   +  B3^3   +  B4eS4   +  B5e3n5   +  B^6    +  B^r,7 
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If  K2n  is  the  space  of  complex  polynomials  of  degree  n  and 
D(k,n  -  k)  is  the  manifold  of  polynomials  R2n  having  k  roots  to 
the  left  and  n  -  k  roots  to  the  right  of  the  imaginary  axis  of  the 
complex  sphere,  then  "by  setting  up  the  following  table: 


AQ  A±     A2  A3  A4  A5  A6  A?  Af 


0 


Bl     B2     B3     B 


B5     B6     B7     B6 


(3.6) 


and  by  making  the  transformation 

%  +  \lBl       A1  +  \lB2       A2   +  \b3        A3   +  \lB4. .  .A~ 
0  Bl  B2  B^      ...C 


(3.7) 


table  (3.7)  is  found  to  correspond  to  a  polynomial  of  the  same  type 
with  respect  to  the  distribution  of  the  roots  relative  to  the  imaginary 
axis,  as  in  equation  (3.4). 

From  table  (3.6),  an  inequality  is  obtained  that  defines  the 
region  in  the  plane  KN  corresponding  to  the  presence  of  the  first 
root  of  equation  (3.1)  in  the  upper  half -plane: 

AqB-l  <  0  (3.8) 

By  setting     \^  =  -  Aq/B-j_     in  table    (3.7) 

(A1B1   "  Vz^l    (A2E1   "  A0B3}/B1    (A3B1   "  A0B4)/B1    —  \  A8"l 


B-, 


B, 


B, 


..B7   0 


(3.9) 


Because     A1B]_   -  A^B2   =  "   16K(K  -   l)3  <    0     for     K>1,   by  multi- 
plying the  elements  of  the  first  rows  of    (3.9)  by     B-L2/(A1B1  -  ApBg) 
and  changing  signs  in  the  second  row 


B- 


-    Bn 


Dl 

B2 


-  B3      -  B4 


D„ 


Bc 


Dc 


D£ 


B6      -  B7 


(3.10) 
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wnere 


Bl<Vl   -  A0B5) 
A1B1   "  ^2 

B1(A5B1   -  ACtBA) 
A1E1   "  V2 


=  Dl 


=  D2 


B1(A4B1   -  A$5)    _ 
A1B1   "  AcP2 

Bl(A5Bl   "  A0B6) 


=  D, 


A1B1  - 

^2 

Wi 

"  AoV 

A1B1  - 

A0B2 

B-jA 

7 

A1B1   - 

A0B2 

B±A 

8 

=  Da 


=   He 


(3.11) 


A1B1  "  A0B2 


The  first  row  of  table  (3.10)  is  left  unchanged  tut  to  the 
second  row  is  added  the  first  row.  Thus 


Bn 


Dn 


D. 


D, 


D, 


Dc 


D£ 


D71 


0   Dx  -  B2  D2  -  B3  D3  -  B4  D4  -  B5  D5  -  Bg  Dg  -  B?   C? 

(3.12) 

From  the  preceding  calculations,  an  inequality  is  obtained  that 
defines  the  region  in  the  plane  KN  that  corresponds  to  the  presence 
of  the  second  root  of  equation  (3.1)  in  the  upper  half -plane 


B1(D1  -  B2)<0 


(3.13) 
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By  carrying  out  a  transformation,  similar  to  (3.7)  of  table  (3.12), 

B-l  +  ^2(D1  "  B2^   Dl  +  ^2(D2  "  B3^   D2  +  ^2^3  ~  B4)...D^ 
0  D,  -B0  D0  -  B^     ...C- 


Dl  "B2 


D2  "B3 


'7 


K3.14) 


By  setting  Xg  =  ~  B-j_/(D-,  -  B2)   and  substituting  in  (3.14) 


D1(D1   -  B2)    -  B1(D2   -  B3)        D2(D1   -  Bg)    -  B1(D3   -  B4) 


Bl   "  B2 
Dl   "  B2 


%   -  B2 
D2    -  B3 


^(Di   -  B2)    -  B!(D4    -  B5) 


>    (3.15) 


Bl   "  B2 
D3   -  B4 


Because     D1(D1   -  B2)    -  Bx      (D2   -  B,)  >  0     for     K>1,   by  multi- 
plying the  elements   of  the  first  row  of    (3.15)   by   (D-j_   -  B2)    / 
[D^Di   -  B2)    -  B^Dg   -  B3)] 


( 


%   ~  B2 


[P2(D1   -  B2)    -  B1(D3   -  B4)]  (Dx   -  B2) 
DiCD-L   -  B2)    -  B1(D2   -  B3) 


LD1   "  B2 


D2   -  B3 


.] 


(3.16) 


The  elements  of  the  first  row  are  subtracted  from  the  elements 
of  the  second  row  of  table  (3.16). 


D-l  -  B2 


[(D2  (Dx  -  B2)  -  B1(D3  -  B4)]  (D1   -  B2) 


0    (D2  -  B3) 


[DgCDi  -  B2)  -  B1(D3  -  B4)]  (Dx  -  B2) 


1 


D1(D1  -  B2)  -  B1(D2  -  B3)     "J 

(3.17) 

From  the  preceding  table,  an  inequality  is  obtained  that  defines 
the  region  in  the  plane  of  the  parameters  KN  that  corresponds  to  the 
presence  of  the  third  root  of  equation  (3.1)  in  the  upper  half -plane. 


(Di    -  Bo) 


(Dg   -  B3)    - 


[D2(DX   -  B2)    -  B!(D3    -  B4)]  (Di    -  B2) 


D1(D1   -  B2)    -  B1(D2    -  B3) 


<0 
(3.18) 
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Similar  conditions  can  "be  obtained  for  all  the  remaining  roots 
of  equation  (3.1).   This  investigation  has  been  limited  to  the  three 
conditions  that  are  sufficient  for  proving  the  existence  of  several 
unstable  waves. 

By  replacing  inequalities  (3.8),  (3.13),  and  (3.18)  by  equations, 
the  equations  of  the  curves  determining  the  breakup  of  the  KN  plane 
into  regions  are  obtained.  The  most  interesting  case  of  large 
K  =  n-)_/(i2^>l  is  considered.   From  inequalities  (3.8),  (3.13), 

and  (3.18)  and  by  considering  equations  (3.2)  and  (3.11)  and  neg- 
lecting small  powers  of  K,  the  following  equations  are  obtained: 

2(K  -  l)3(K  +  3)  =  0 
egW3  +  e-]_N2  +  e2N  +  e3  =  0 

4(K2  -  1)(K  +  3)N  +  K(K  -  l)(K3  +  17K2  -  96K  +  99)  =  0 

(3.19) 


where 


el 


e0  =  128  (K4  -  K3  -  23K2  -  39K  -  18) 
=  592K(K5  +  8.4K4  +  3.18K3  -  96K2  -  20. 3K  +  0.98) 


e2  =  9K2(K6  +  8.4K5  -  97. 3K4  -  2045K3  +  1700K2  +  390K  +  363) 
e3  =  24K5(K5  +  12. 3K4  +  306K3  -  4100K2  +  12,300K  -  7000) 

By  plotting  the  curves  (3.19)  in  the  KN  plane  and  separating 
by  hatched  lines  the  regions  corresponding  to  the  signs  of  the 
inequalities  (3.8),  (3.13),  and  (3.18),  the  diagram  shown  in 
figure  3  is  obtained.  This  diagram  shows  that  for  K>  0  and  N>0 
a  region  of  values  of  K  and  N  exists  that  corresponds  to  the 
presence  of  three  roots  with  positive  imaginary  part,  that  is, 
of  three  unstable  waves  on  the  surface  of  separation. 

The  division  of  the  KN  plane  for  the  remaining  roots  could 
establish  regions  with  a  still  greater  number  of  roots  with  posi- 
tive imaginary  part.  The  given  incomplete  diagram  already  shows, 
however,  the  existence  of  several  unstable  waves.  In  the  presence 
of  a  maximum  B^  or  c^,  several  infinitely  long  strings  will 
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break  away  from  the  surface  of  separation,  the  cross-sectional 
dimensions  of  which  will  depend  on  the  wavelength  of  the  critical 
disturbance. 


Translated  by  S.  Reiss, 
National  Advisory  Committee 
for  Aeronautics. 


REFERENCES 

1.  Rayleigh:  The  Theory  of  Sound.  Dover  Pub.,  2d  ed.,  1S45. 

2.  Tomotika,  S.:  On  the  Instability  of  a  Cylindrical  Thread  of  a 

Viscous  Liquid  Surrounded  by  Another  Viscous  Fluid.   Proc. 
Roy.  Soc.  London,  vol.  CL,  no.  A870,  ser.  A,  June  1935, 
pp.  322-337. 

3.  Petrov,  G.  I.:   On  the  Stability  of  Turbulent  Layers.   Rep. 

Wo.  304,  CART,  1937. 

4.  Neimark,  Y.  I.:  On  the  Problem  of  the  Distribution  of  the  Roots 

of  Polynomials.  DAN,  T.  58,  No.  3,  1947. 


NACA  TM  1281 


19 


vt.Pi.N 


mm. 


WW////M 


Figure  1. 


Figure  2. 


NACA-Langley  -  4-16-51  -  875 


(O 

to 

C 

• 

o 

H 

•H 

• 

-P 

<-\ 

NACA^ 

es  on  Surface  of  Separa 

•H 

-P 
•H 

n 

• 

/     > 

• 

>H 

r      & 

CO 

f       > 

T3 

-O 

>j 

•H 

3 

Cl 

u 

rH 

<S 

P-H 

Pi 

H 

•H 

CO 

H 

CO 

■g 

fl 

•H 

3 

•H 

Pi 

o 

k 

rH 

M 

aJ 

o 

o 

CD 

•H 

o 

co 

w 

00 

rH 

2 

•M 

rH 

ir\ 

03 

> 

• 

ON 

-p 

.-H 

< 

S 

r-\ 

<D 

r^ 

o 

EH 

•"3 

p 

£ 

> 

< 

r-i 
•rH 

•\ 

tn 

o 

Jh 

g 

O 

pq 

gj 

ft 
< 

rH 

Ph 

Cm 

<Q 

O              >5    to     i>J 

^ 

1    ro             1-1    ©  -d 

•P 

5-,  -h    ©         ©    f>    3 

O 

r* 

d 

3         0        -p    ro  -P 

>    © 

CO 

•M 

ffl-P    C        -h   £    ro 

-P  -P 

0 

-P 

3 

Pi   J)     O    Cl 

ro 

c 

C* 

«    S    inflH    !>iO 

Cm    5-i 

ro 

Cm 

© 

X    -M 

Si    ©    CO"   -P   *H    Jh   40 

O  -P 

© 

O 

M 

tSDrH 

Ch 

-P  .p    ©           C    (S 

to 

> 

© 

=! 

O 

43     ft(H    -H   H     © 

M      fl 

?3 

© 

Ch 

0  ro 

C     (J      PlO             rH    r-l 

O      O 

^5 

O 

Cm 

M 

ro 

O             CS             d    "H    ,£>. 

-h  a 

ro 

•H 

cm 

ro 

0,            .p    -M     Pj  -M 

+3  © 

© 

Ch 

•d 

ro  0 

© 

CO    03     ©   rH           ro     to 

ro  <d 

ro 

H 

0 

©          ^    S-P    O    B 

5-1 

© 

3 

ch 

ro  a 

O 

f>      «*-p    to    ©           O 

ro  0 

X 

CO 

O 

©      M 

M 

03   ro        ©  tj  ©   a 

ft-P 

FM 

>  © 

P4 

>    -d     rH     rH            ^ 

© 

a 

d 

•H    -p 

•HO             ©    +3     CD 

ro  >d 

c 

■H 

to-p 

© 

>1JIH     Jlfl            «H 

H 

• 

■H 

3 

ro 

X 

5n    cr*         ,d+"H 

Cm   ro 

ro 

-P 

H 

>d   P,  +3 

03   «H     to  +3           O   -P 

0 

© 

•~H 

Cm 

© 

rH    r-H    T-t            +3            -H 

ro 

> 

-P 

-P     CJ 

c! 

rH               t0     CQ     03      CD 

©  -d 

03 

u 

Cm 

ro  0 

•H 

■ri     51     OJ     (S^^       *> 

O    -H 

JS 

ro 

O 

Cfl-H 

d 

-p 

Pi     3     rQ               +3     -P     -P 

ro  3 

p< 

•H    P> 

H 

o 

c3    O          to          CO  © 

Cm    rH 

>s 

CO 

-p  ro 

a 

cri 

O     O    rH      ?n     t0    £    t-J 

m  Cm 

U 

© 

M 

ro   5-1 
©  ho 

R 

U 

ro    ro    ©    C    © 

3 

ro 

,3 

s 

<D 

■P 

©    iH     O    -P    -H    rH    d 

to    to 

rH 

•p 

•H 

t>   © 

0: 

rH     t>    iH     ©     S            -H 

3 

rH 

u 

C  P 

O 

rO 

rQ            -P     33     ©     3 

©    O 

•H 

E 

■p 

•h  a 

-P 

< 

c3    0    03    ro    to  ^   o1 

C    0 

P< 

O 

m 

•M 

+3   ts   a  -h   to  -p  -h 

ro  to 

cd 

fn 

S   ro 

d 

tO    -P      ©    T3      ro             rH 

r-{    -H 

O 

«H 

tfi 

©    -H 

a 

a      si          0 

P<    f> 

C 

rH    rd 

© 

S^-PP    1>j-P    © 

83 

>J 

0 

rQ 

-P 

0  ro   a  pq       ,3 

©     C£)rH 

ro 

rH 

O      © 

© 

to        S    ©         C  P 

x!  a 

X 

> 

r(    X 

M 

M      M                5-1               O 

+3    -H 

c3 

§ 

!>a 

Pl-P 

P. 

•H    O    03    ©      •    CQ   Cm 

-d 

+J 

x 

rH 

5-1   "H            <+H    -P    "H     O 

Cm    c! 

to 

cS 

ID 

©  cm 

-P 

©    +3       ©     Ch      O      ^ 

O    © 

ri 

© 

-P 

Si      O 

O 

•d  ro  d  -h  '-j  ro  © 

>>"£! 

3 

rl 

■H 

Eh 

a 

■H     r<    "H    -d              Pi    O 

rQ 

a. 

Cl 

ro  ro  >       ro  a  ro 

-P    © 

Cm 

■H 

O 

ro 

• 

fl    a    O  Cm           0  Cm 

•H 

O 

© 

Cm 

•    -rl 

© 

a 

O    ©     M     O   cm     O     M 

rH     !>5 

S\ 

a 

CO   -p 

0 

0 

O    to    PL,         0          3 

•M    rH 

© 

-P 

•H 

m  ro 

d 

•M 

to         era 

rO      © 

O 

0  a 

-P 

S^Ch    0  -P    Pj  -h 

ro  -P 

S 

O 

rH 

•M    -M 

+3 

ro 

pq    O   +3    ©    3          © 

-p    -H 

© 

-P 

03 

to    M 

3 

N 

rH  ,M    ©  ,d 

co   S  -p 

M 

c  0 

rQ 

■M 

©  ©  ft  ro  to-p 

•H 

m 

>d 

© 

©     M 

a 

0  d    O    ©    U 

rM^ 

•H 

ro 

> 

a  p.-p 

0 

ro  ro  m  ?h  ro  £ 

H 

© 

© 

•H    p. 

© 

-p 

Cm    0  <d  rQ  rH    O 

-P   -H 

© 

rH 

ro 

-d  ro 

T-3 

ro 

UNIVERSITY  OF  FLORIDA 


3  1262  08105  058  4 


'  OF  FLORIDA 

"MENT 

;CE  LIBRARY 
17011 

VILLE.FL  32611-7011  USA 


